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ABSTRACT 

We present a model of an outburst of the soft X-ray transient A0620-003. A two- 
dimensional time-dependent smoothed particle hydrodynamics scheme is used to sim- 
ulate the evolution of the accretion disc through a complete outburst. The scheme 
includes the full tidal potential of the binary and a simple treatment of the thermal- 
viscous disc instability. In the case where the mass accretion rate onto the primary 
determines the fraction of the disc that can be kept in a hot, high viscosity state by 
the resulting X-ray emission, we find that the shape of the X-ray light curve is ulti- 
mately determined by the relative sizes of the irradiated and unirradiated parts of the 
disc and the growth time-scale of the tidal instability. The model accounts for the re- 
brightening that has been observed in the light curves of A0620-003 and several other 
transients. The primary maximum and subsequent decline are due to the accretion of 
gas within the irradiated portion of the disc, while the secondary maximum is caused 
by the accretion of gas in the outer part of the disc that is initially shadowed from the 
central X-rays, but subject to tidal forces. We propose that tidal effects at the disc 
edge can be sufficient to drive accretion on a time-scale shorter than that expected 
for a standard alpha-viscosity disc. The final decay is subsequently controlled by the 
gradual retreat of the irradiated portion of the disc. If the entire disc is kept in the 
high-viscosity state by the irradiation, no rebrightening is possible. 

Key words: accretion, accretion discs - instabilities - binaries: close - X-rays: stars 
- methods: numerical - stars: individual: A0620-003. 



1 INTRODUCTION 

The soft X-ray transients (SXTs) are close interacting bi- 
nary systems consisting of a black hole or neutron star pri- 
mary and a low-mass Roche-lobe-filling secondary star. The 
outbursts of soft X-ray transients are quite distinct from 
those seen in other binary systems such as dwarf novae. The 
duration and recurrence times are both strikingly longer, 
suggesting that a large fraction of the mass of the accre- 
tion disc is consumed in each outburst and that a different 
mechanism from the standard disc instability picture is re- 
quired (comprehensive reviews of the disc instability picture 
are given by Cannizzo, 1993 and Lasota, 2001). 

Typical ionization instability-driven dwarf nova out- 
bursts last a few days and consume a few percent of the disc 
mass, leading to a recurrence time of the order of weeks. 
In comparison, a typical SXT outburst lasts for 100 days or 
more and recurs on a time-scale of several years. The profiles 
of the X-ray hght curves of SXT outbursts vary greatly but 
they do share common characteristics - the initial peak is 



E-mail : mtr@star.le.ac.uk 



often followed by an exponential decline and, in many cases, 
a rebrightening several tens of days after the first maximum. 
Smaller rebrightenings later on in the decay from outburst 
have also been observed. These features can be seen quite 
clearly in the X-ray light curve of the well-known transient 
A0620-003 (Elvis et al. 1975, Whelan et al. 1977). A general 
review of SXT outbursts is given by Tanaka & Lewin (1996). 

Several authors have proposed models to explain the 
form of the SXT light curves. The long recurrence time of 
the outbursts is a problem which has been addressed in many 
ways. It can be explained by the standard ionization insta- 
bility model if the Shakura-Sunyaev viscosity parameter a 
is greatly reduced from that typical of a cataclysmic vari- 
able disc. Such an explanation is not confined to SXTs; it 
has been invoked to explain the outbursts of WZ Sagittae 
(Smak, 1993). However, there is no obvious physical reason 
why a should vary so much from one accretion disc to an- 
other. The more successful models involve the removal of the 
inner part of the disc where the viscous time-scale is short. 
This can been achieved by a hot, advection-dominated accre- 
tion flow (ADAF) onto the hole : Meyer & Meyer-Hofmeister 
(1994) showed that gas from the cool outer accretion disc in 
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a binary system may evaporate into a hot corona and sub- 
sequently feed an ADAF onto the comapct object. Dubus 
et al. (2001) have performed one-dimensional simulations 
of an SXT including the effects of both irradiation and an 
ADAF. In this work wc do not address the question of recur- 
rence time, but concentrate on the observed outburst pro- 
files. Cannizzo, Chen & Livio (1995) found that the expo- 
nential decay could be reproduced by parameterising a in 
terms of the aspect ratio S . Augusteijn, Kuulkers & Shaham 
(1993) sought to explain the rcbrightenings by considering 
irradiation of the secondary as a mechanism for increasing 
the mass transfer rate. However, it is unlikely that the ra- 
diation can penetrate the photosphere of the companion to 
produce a variation in mass transfer rate of more than a few 
percent (King, 1989). 

We consider the effect of disc irradiation by the X-ray 
emission from the compact object. King & Ritter (1998) 
showed analytically that the assumption that the central 
X-ray source can efficiently irradiate the disc leads auto- 
matically to long exponential X-ray decays with the ob- 
served time-scales. Observations of SXTs suggest that the 
optical/X-ray flux ratio for these objects is simply too large 
for the optical emission to result from intrinsic viscous dissi- 
pation in the disc alone (van Paradijs & McClintock, 1994; 
Shahbaz & Kuulkers, 1998). Moreover, the optical brightness 
is strongly correlated with the irradiating X-ray flux. How- 
ever, detailed one-dimensional modelling of the disc struc- 
ture (e.g. Cannizzo 1994) suggests that the response of the 
disc to the irradiation is for the surface to become convex 
rather than concave, thus shielding the outer regions from 
the irradiating flux. The observations seem to contradict this 
picture, and for the purposes of this paper we assume that 
the disc surface remains concave. 

In section 2 we describe the irradiated disc model of 
King and Ritter in more detail, and outline the importance 
of tidal forces in systems with extreme mass ratios in section 
3. The rmmerical method used to implement the model is 
then described in section 4. In section 5, wo show that the 
observed peak luminosity of the 1975 outburst of A0620- 
003 almost certainly implies that the irradiation radius (the 
outer radius of the region kept heated by the central X- 
rays) was larger than the accretion disc. We then show in 
simulations that no rebrightcuiug occurs if the whole disc 
is kept in the hot viscous state by the irradiation; however, 
if only a very small area of the outer disc is shielded in 
some way (by a radiation-induced warp or self-shadowing, 
for example), the observed rebrightening is produced. 



2 DISC IRRADIATION MODEL 

Irradiation of an accretion disc by the central X-ray source 
heats the exposed disc faces. In SXTs, provided the disc 
faces remain concave, the surface temperature of the disc is 
given by the formula of de Jong, van Paradijs & Augusteijn 
(1996) : 
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where rj is the accretion efficiency for producing X-rays, M 
is the mass accretion rate onto the primary, H is the scale- 
height of the disc and /? is the disc albedo. The index n 
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Figure 1. Schematic representation of a SXT outburst according 
to the King & Ritter (1998) model. 



takes the values 1 or 2 for irradiation by a point source or 
by the inner disc respectively. Accretion onto a neutron star 
will always have n = 1, while n = 2 is possible in the ab- 
sence of a hard surface on which to accrete, i.e. a black hole 
primary. During outburst, however, the disc will probably 
develop an X-ray corona, as suggested by observations of 
a strong power-law continuum. In this case, the appropri- 
ate index reverts to n = 1. In a large enough disc, Tirr is 
quite sufficient to dominate the effective temperature and 
exceed the hydrogen ionization temperature (Th ~ 6500K), 
thus stabilising the disc on the hot branch of the Tefi — E 
relation (sec for example, Cannizzo, 1993). Tin is the key 
factor which determines whether a low mass X-ray binary 
is transient or not - in a persistent source Tin exceeds Th 
throughout the entire disc (van Paradijs, 1996). 

We consider a situation in which the disc undergoes an 
outburst which is initiated by the standard hydrogen ion- 
ization disc instability. Subsequently, as the mass accretion 
rate onto the primary increases, the disc is irradiated by 
the enhanced X-ray emission. The disc is unable to return 
to a cool quiescent state until the mass accretion rate has 
fallen sufficiently to switch off the enhanced irradiation, al- 
lowing a cooling front to sweep through the disc. Following 
the analytic treatment of King & Ritter (1998), one can use 
equation 1 to determine the variation of Rh with Mc by not- 
ing that the edge of the irradiated region is determined by 
the point in the disc at which Tin = Th. Then : 
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Using typical parameters for neutron star and black hole 
systems, King & Ritter showed that the mass accretion rates 

for a given Rh arc then 

IOt 



Mi.crit = 4.1 X 10"'"Rii Moyr" 



(3) 



and 

M2,crit = 2.9 X 10-^Rii^Moyr-\ (4) 

where Rn is Rh in units of 10^^ cm. King and Ritter and 
subsequently King (1998) extend this analysis to show that 
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the first part of tiie outburst decay, witfi tiie wliole disc irra- 
diated, wiii be exponentiai, but uitimateiy tlie finai part of 
the decay will be roughly linear due to the gradual shrink- 
ing of the irradiated region. However, they could not fully 
explain the rebrightening, and its origin remained unknown. 
A schematic of the theoretical profile is given in fig. |^. 

The recurrent transient A0620-003 (V616 Mon) is per- 
haps the best known of the SXTs. Its 1975 outburst was 
observed by Ariel V (Elvis et al. 1975), during which sec- 
ond and third maxima were observed 55 and 200 days after 
the first. The decay either side of the first rebrightening is 
exponential and the rise-time of the first rebrightening was 
~ 10 — 15 days. A0620-003 has a short orbital period of 
7.75 hours (McClintock et al., 1983), and is likely to contain 
a black hole primary. The primary mass has been estimated 
to be 7.3M0 or more (McClintock & Remillard, 1986), but 
a more recent work concludes that it is likely to be in the 
range 3.30 < Mi < 4.24 M© (Marsh, Robinson & Wood 
1994). The reader is referred to Whelan et al. (1977) or 
Tanaka & Lewin (1995) for the full X-ray light curve of the 
1975 outburst. 

The analysis of King & Ritter correlates well with the 
observed behaviour of A0620-003 and SXTs in general, but 
the precise origin of the rebrightening and the behaviour 
of the disc remain to be explained. In this work we incor- 
porate King & Ritter's analytic model into an existing nu- 
merical code and perform two-dimensional simulations of a 
SXT with the orbital parameters of A0620-003. We find a 
natural explanation for the rebrightening in terms of the 
physical state of the disc. In a forthcoming paper we will 
report results from simulations of outbursts in longer period 
SXTs (with correspondingly larger discs) and investigate the 
resulting complex outburst morphologies. 



4 NUMERICAL METHOD 

We use a two dimensional smoothed particle hydrodynam- 
ics (SPH) accretion disc code. The great advantage of this 
method is that it incorporates the full tidal potential of the 
binary system without need for an approximate prescription 
for the removal of angular momentum at the disc edge em- 
ployed in 1-D simulations. The importance of using a two- 
or three-dimensional scheme with the tidal effects included 
has been demonstrated effectively with this code in recent 
simulations of cataclysmic variables (Truss et al. 2001). Full 
details of the code are described in Truss et al. (2000) and a 
complete review of the SPH technique is given by Monaghan 
(1992). The implementation here is very closely related to 
the dwarf nova study of Truss et al. (2000,2001), in which 
we incorporate the disc instability in terms of a viscosity 
switch based on surface density and build up the disc by a 
mass transfer stream. The disadvantage of the scheme is that 
there is no self-consistent calculation of vertically-averaged 
thermal structure, as in one-dimensional codes. There is a 
trade-off in computational efficiency between calculating the 
thermodynamics in a self-consistent manner in 1-D and per- 
forming the hydrodynamics in 2-D. We choose to present 
a 2-D calculation in this work to examine the role of tidal 
forces, and in this way the scheme is complementary to mod- 
ern 1-D treatments. 



4.1 Viscosity 

When the local surface density at the position of any particle 
in the disc exceeds the critical value for ionization of hydro- 
gen to occur as calculated by Cannizzo, Shatter & Wheeler 
(1988), 
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3 TIDES 

The mass ratio of A0620-003 is very well-constrained : q = 
^ = 0.067±0.01 (Marsh. Robinson & Wood, 1994). Such a 
low value suggests that tidal resonances will be accessible to 
the disc. Tidal forces will therefore be an important factor 
in the dynamics of the accretion disc. The radius of the 3:1 
tidal resonance, which has been shown to be important in 
dwarf novae, is 
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where a is the binary separation (see, for example, Frank, 
King & Raine, 2002). For q = 0.07, this radius is R3:i = 
0.47a. This is well inside the radius at which the disc is 
ultimately truncated by the tides. Rtidc can be estimated 
from the Roche lobe radius of the primary - it is always 
roughly 90% of this value. Using the formula of Eggleton 
(1983), 
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Rtide — 0.60a for q = 0.07. Therefore we expect tidal effects 
to become important and the disc to become eccentric when 
it reaches Rsii. 



the viscosity acting on that particle is increased. The change 
in viscosity is made on an appropriate thermal time-scale 
(much shorter than the local viscous time-scale) from a value 
Qcoid to a value ahot, where a is defined in terms of the 
sound speed of the fluid, Cg, and the scale height of the disc, 
H, in the usual viscosity formulation of Shakura & Sunyaev 
(1973): 

V = |aCeH. (8) 

The local viscosity remains at dhot until the local surface 
density decreases below a critical value 
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when it is returned to the low state value QcoM- 

Equation ^ is employed almost directly in the SPH for- 
mulism. The SPH momentum equation for a particle i is: 



dvi >r-^ /Pi Pj Pfifj - CcijAtij 
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where pij and Cij are the average density and sound speed of 
the pair of particles i and j respectively, h is the smoothing 
length and Wij is the interpolating kernel. /3 and (" are the 
quadratic and linear artificial viscosity parameters and 
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Figure 2. Radial variation of SPH viscosity and a-viscosity. The 
two solid lines show the viscosity in an a-disc with a = 0.1 and 
Ml = 4Mq. The upper is for M = IO-^Mqjv-'^ and the lower 
for M = lO~^M0yr~^. The dashed lines show the SPH viscosity 
where C = 1, Mi = 4M0, Cs = O-OSaHb and C = 1-5 (top), ( = 1.0 
(middle) and = 0.5 (bottom), a is the binary separation and Slj, 
is the angular velocity of the binary. 



where 77^ = O.Olh'^ avoids singularities. 

Following Murray (1996), we set j3 — and in the con- 
tinuum limit (i.e. h — > and number of particles N — + 00) 
the acceleration due to the artificial viscosity term can be 
written 
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where k is a constant determined by the kernel W. In this 
case we use a 2D cubic spline kernel, for which «: = i. Sij is 
the deformation tensor, where 



Si' = — + ^ 



(13) 



If the density and sound speed vary on lengthscales much 
longer than the velocity, and if we assume that the fluid is 
divergence-free, then 



(14) 



So, in the continuum limit, the artificial viscosity introduces 
a shear viscosity 

u = iCcsH (15) 

In fact, Murray (1996) used the SPH smoothing length, h, 
as the lengthscale in his prescription for viscosity. We do 
not impose such a constraint here : the viscous energy is 
dissipated over the disc scale height, H. In the simulations, 
we take the sound speed to be constant everywhere and vary 
1^ to change the viscosity of a particle in the way described 
above. With such an isothermal scheme, the net viscosity in 
the simulations has a slightly different radial dependence to 
an analytic alpha-disc. In the alpha model, the variation of 
sound speed and scale height with radius (cs ^ R^s and 
H ~ rI) lead to 



where Mie is the mass accretion rate in units of 10^ gs~^. 
In the SPH simualtions, the sound speed is a constant and 
the scale height is a linear function of radius: 
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where we usually use ^coid = 0.1 and C^hot = 1-0. A com- 
parison of these two viscosity variations is shown in fig. ^ 
A full discussion of how the results of this scheme can be 
interpreted in relation to observations, and its limitations, 
is given in Section 6. 

4.2 Radiation 

The irradiation radius is calculated by direct implementa- 
tion of equation ^ at each timestep . All particles within 
this radius are instantaneously given the high-state SPH 
viscosity (^hot. Since we know the number of particles being 
accreted onto the primary in the simulation, we can readily 
calculate the mass accretion rate at each timestep, and 
hence the irradiation radius. We are not concerned with 
calculating accurate optical light curves in these simula- 
tions; instead we concentrate on the total X-ray emission, 
which we take to be proportional to the mass accretion rate 
onto the primary. This scheme is ideal for investigating how 
variations in the X-rays arise from waves of mass travelling 
through the disc. 

Our model is distinct from those employed by Huang 
& Wheeler (1989) and Mineshige & Wheeler (1989) in that 
we follow the evolution of the whole disc in two dimensions 
and there is no need to approximate the binary potential. In 
particular, for numerical reasons, Mineshige & Wheeler used 
a very small disc in their simulation of A0620-003, using an 
outer radius of 3.16 xlO^^cm, well inside even the circulari- 
sation radius. We simulate the entire disc and show that the 
outermost regions can play a pivotal role in determining the 
form of the observed X-ray light curve. 



5 RESULTS 

5.1 Irradiation efficiency 

Given a distance of 870 pc (Oke, 1977), the peak X-ray 
luminosity of the 1975 outburst of AQ620-003 was Lx ~ 
1. 3x10"^* erg. s~^ . The total accretion luminosity can be writ- 
ten 



: 77MC 



(18) 



where 77 is the rest mass to energy accretion efficiency. We 
can use this equation with King & Ritter's analytic expres- 
sions for Rh to estimate the maximum value of Rh reached 
at the peak of the 1975 outburst. Equations (3) and (4) were 
derived using specific values of rj; (3) assumes that the effi- 
ciency for accretion onto a neutron star is 77 = 0.15, while 
for the black hole case (4) assumes ri — 0.10. If we relax this 
constraint and rewrite the equations in terms of a general 77 
consistent with ([18|), we have : 



Rf 



2 -1 
cm s 



(16) 



Mc 



and 



6.2x10" 
V 



-R?i Moyr"' 



(n=l), (19) 



SXT outbursts 5 



3.5 r 




n.n h I , . . .1 I I. . . I I I i 

2 3 4 5 6 7 8 

M,+M, (M,„|) 



Figure 3. The upper curve is tlie binary separation, a, for 
Ml = 3.7M0 and the lower curve is the Roche lobe radius of 
the primary, Rl(1), in A0620-003 with a mass ratio q = 0.07. 
The disc radius will not exceed this value. The range of masses 
derived by Marsh, Robinson & Wood (1994) is denoted by vertical 
lines and the dashed horizontal line shows the maximum value of 
Rli deduced from the observed peak luminosity of the 1975 out- 
burst and the weaker (n=2) prescription. It is clear that Rjj is well 
outside the disc and that this constraint does not change even for 
masses in excess of 7M0. The n=l prescription leads to a radius 
of Rji = 6.1xl0^^cm, well above the upper limit of the plot. 

Merit = R?i Meyr-i (n = 2). (20) 

n 

Now, if we infer the mass accretion rate from the ob- 
served luminosity in ([l8|), we see from ( [l9| ) and (^o|) that 
the estimates of Rh are independent of r). For n=l we infer 
that Rh = 6.1 X 10"cm and for n=2, Rh = 2.8 x 10"cm. 
For a total system mass of Mi -I- M2 = 4M0 and orbital 
period 7.75 hours, the binary separation of the system is 
2.2 X 10^^ cm. Therefore, the maximum value reached by Rh 
during an outburst of A0620-003 is outside the disc, and 
quite possibly much larger than the disc. The parameter 
space for varying masses is shown in figure ^ Given a mass 
ratio, q, the disc size can be estimated using the formula for 
the volume radius of the Roche lobe of the primary, equa- 
tion ^. This is shown in the lower curve of fig. |^, for q = 
0.07. 



5.2 Outbursts 

We perform separate simulations of A0620-003, allowing the 
irradiation radius to take different values for the same mass 
accretion rate. In this way we are effectively varying the 
parameter Bn in equations. We aim to discover the effects of 
allowing Rh to become larger than the accretion disc or not. 
In these simulations we take Pd-b = 7.75h, Mi = 3.7M0 and 
M2 = O.25M0, consistent with the measurements of Marsh, 
Robinson & Wood (1994) and -M2 = 3xlO""M0yr-^ after 
McClintock et al. (1983) and Cannizzo (1998). The sound 
speed is Cs — 0.05af2b ~ 24.5kms~^ and the high and low- 
state SPH viscosities are (^hot = 1 and C^coid =0.1. 

Figure |^ shows the variations in mass accretion rate 
onto the primary and heating radius Rh for the simulation 
of A0620-003. The entire disc is irradiated and there is no 
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Figure 4. Outburst simulation of A0620-003 with Rh,max >> 
Rdisc- The upper panel shows the mass accretion rate onto the 
primary. Below is the variation in Rh through the simulation. The 
decay is initially exponential, but when the irradiation retreats 
inside the disc, the decay becomes faster. 

large rebrightening. There is an exponential decay as the 
irradiated gas is accreted followed by a faster decline as the 
accretion rate is controlled by Rh gradually retreating inside 
the disc. There is a small variation at 70 days after outburst 
maximum. This arises from the hot / cold boundary passing 
through structure near the disc edge associated with the 
spiral arms. The effect is too small to be responsible for the 
larger rebrightening in A0620-003, but may be responsible 
for the later, smaller glitches observed in the light curve. 

In the second simulation, shown in fig. ^ the efficiency, 
ri, is reduced such that Rh is limited to a maximum value 
Rh.max = 0.35Rdisc- This allows us to mimic the effect of 
a self-shadowed outer disc. In this case, when the mass ac- 
cretion rate starts to rise, only the inner part of the disc 
is irradiated. The accretion of gas in the irradiated region 
gives rise to the first peak in fig. ^, which decays exponen- 
tially on a time-scale appropriate to the inner part of the 
disc. Meanwhile, the disc encounters the 3:1 tidal resonance 
and a two-armed spiral wave is launched in the disc. Gas in 
the outer, unirradiated regions is accumulated in the spiral 
arms and the surface density rises. Eventually, E exceeds the 
critical limit for the viscosity to be increased. When this gas 
subsequently arrives at the primary, there is a rebrightening 
in the X-ray light curve. The second stage of the decay takes 
place on a longer time-scale, and the slope is correspondingly 
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Figure 5. Simulated outburst of A0620-003 with t] reduced so 
that Rh.max = 0.35Rdisc- The point at which Rh becomes smaller 
than 0.35Rdisc is marked with an arrow. 
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Figure 6. Run A : Simulated outburst of A0620-003 with rj 
reduced so that Rh,max = 0.85Rdisc- The arrow marks the point 
at which Rj^ becomes smaller than 0.85Rdiac- 

more gradual. The final, rapid drop to quiescence only takes 
place wlien Rh < 0.35Rdisc- 

Tiie two distinct regions, one irradiated and one unir- 
radiated but subject to tiie disc instability, give rise to the 
rebrightening. The light curve looks more like that of A0620- 
003, but differs in one important detail - the slope of the 
decay either side of the rebrightening is not the same. The 
observed similarity in the gradients of the decay either side 
of the rebrightening in A0620-003 suggests that the waves of 
mass which are producing these parts of the outburst origi- 
nate in the same region of the disc. 

It is clear from the discussion in section 4.1 that Rh 
becomes much larger than the outer edge of the disc dur- 
ing an outburst of A0620-003. However, if a small part of 
the outer disc is shadowed from the radiation in some way, 
the rebrightening and the decay slopes could be reproduced 
even when Rh >> Rdiac- In order to demonstrate this, a 
simulation was performed in which the efficiency limits Rh 
to a value near, but not quite as large as Rdisc (Run A). It 
is not our intention to suggest that the accretion efficiency 
is extremely low, but to perform a simulation that mimics 
the shadowing of a small part of the outer disc during out- 



Run Chot No Figure 



A 


1 


49814 


6 


B 


1 


24996 


7a 


C 


1 


12588 


7b 


D 


1.5 


49814 


8a 


E 


0.5 


49814 


8b 



Table 1. Summary of run parameters. No is the number of par- 
ticles in the disc at the start of the simulation. 
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Figure 7. Outburst simulations with fewer particles. Top : Run 
B. No = 24996 ; Bottom : Run C. No = 12588. 

burst. The resultant light curve is shown in fig. ^, and this 
time the rebrightening and slopes are closer to the observa- 
tions. Again, the final drop to quiescence only occurs when 
Rh < 0.85Rdisc. 

5.3 Varying the simulation parameters 

In order to assess the robustness of these results we per- 
formed run A using different SPH viscosities and different 
particle numbers. A summary of the parameters in each run 
is given in table 1. Firstly, the simulation was repeated with 
the same viscosity as before, ^coM = 0.1, ("hot ~ 1-0, but with 
fewer SPH particles (runs B and C). Figure ^ shows the re- 
sults with roughly one half and one quarter the number of 
particles of the orginal calculation. The greatest difference 
as particle number is decreased is the resolution of the pri- 
mary maximum. With fewer particles, the spatial resolution 
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Figure 8. Top: run D. (^i-^^t = 1-5 ; Bottom: run E. (^i-^^t = 0.5. In 
both runs, Nq = 49814. 



in the inner part of the disc where the density is lower suf- 
fers to a greater extent than the outer part. Consequently, 
the accretion of particles from the inner disc leads to an 
under-resolved maximum in the accretion rate. Therefore, 
we expect the primary maximum to be even better resolved 
for a higher particle number, and the convergence between 
runs A and B is encouraging. The particle masses are scaled 
so that the total mass of the disc is the same in each simula- 
tion, so in run C the large drop in mass transfer rate before 
the rebrightening is caused by the resolution being very poor 
indeed - the accretion of just a few particles corresponds to 
a large response in accretion rate. The accuracy of this run 
should be judged accordingly. 

In runs D and E we keep the particle number the same 
as in run A, but vary the magnitude of (hot- The resulting 
light curves are shown in fig. M. In both runs, the rebright- 
ening starts at the same time as it does in run A - between 
15 and 20 days. However, the rise-time of the rebrightening 
is seen to increse as ("hot is decreased (fig. ^ . This suggests 
very strongly that the the rise-time of the rebrightening is 
governed by the viscosity, but the time delay from the start 
of the outburst to the rebrightening is independent of the 
viscosity. 




Time (days) 

Figure 9. Overlay of the outbursts in run A, with fhot = 1-0 
(solid Hne) and run D, with (^^ot = 1-5 (dot-dashed line). The 
time delay between the onset of the outburst and the start of 
the rebrightening is the same in both simulations; however, the 
rise-time of the rebrightening is longer for the lower viscosity. 



6 DISCUSSION 

In this section we discuss the applicability and limitations of 
these results to interpreting observations of A0620-003 and 
suggest a possible mechanism for the rebrightening. 



6.1 Thermal structure 

We have already seen in section 4.1 that the SPH viscosity 
used in this isothermal scheme has a slightly different func- 
tional form to that of a standard a-disc. Fig. ^ shows that in 
general the viscosity is slightly overestimated at small radii 
and slightly underestimated at larger radii. It should also 
be kept in mind that in the simulations, the irradiated re- 
gion is indistinguishable from an unirradiated region in the 
high-viscosity state. The question thus arises as to how good 
an approximation this is to a real irradiated disc, i.e. does 
the radiation alter the thermal structure? The magnitude of 
the optical depth, r, and the origin of opacity in the disc is 
the key issue. Tirr can be quite sufficient to dominate the 
effective temperature of the disc 



/ 3GMiM(R) 
87raR3 



(21) 



but the vertical thermal structure of an irradiated disc will 
only be altered if Tirr also exceeds the mid-plane temper- 
ature of the disc, given by Tm'' ~ rTeff*. King (1998) has 
argued that t is of order 10 and consequently the vertical 
structure of the outer disc is nearly isothermal, with the 
sound speed dependent on Tirr and not Tm. The viscosity is 
therefore independent of the local value of E and the diffu- 
sion equation 



as 3 9 



dt 



RdR 



leads to 

at 

and the observed exponential decays. 



(22) 



(23) 
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However, simulations by Dubus et al. (1999), in which 
the vertical structure is calculated self-consistently includ- 
ing the effects of irradiation by a point source in the orbital 
plane, suggest that r is of the order 10^ — 10'' and the struc- 
ture is close to that of an a-disc. Indeed, for an a-disc with 
a Kramer's opacity law, 



r = 1200 



m: 



where Mie is the mass transfer rate in units of 10 gs~ 
the more familiar case of an unirradiated accretion disc, 
equation for vertical radiative energy transport is: 



(24) 

i.In 

the 



9 ^GM 



(25) 



and we can see straight away that although the opacity k 
varies only slowly with E and Tm in much of the disc, the 
viscosity i/ is not independent of E. With the alpha prescrip- 
tion, 



9E 
dt 



E3 



(26) 



gives a power-law decay. 

Clearly, the magnitude of r is crucial to determining 
any change in the vertical thermal structure of the disc. It 
remains unclear how the energy is deposited in the vertical 
direction in an accretion disc, and whether the a-model is 
a valid description of the real processes involved. If the a- 
model breaks down we have no way of calculating r. How- 
ever, it remains entirely possible that Tm > Tirr > Th and 
the irradiation ionizes the disc, bringing about a change in 
viscosity, but does not alter the vertical thermal structure. 
The assumption that we make in the code is that the vis- 
cosity increases when Tirr > Th and there is no change in 
Tm (i.e. the sound speed remains constant). 

We do not claim that our approach is a more accurate 
representation than a ID model that is thermodynamically 
self-consistent within the framework of the a-model. Rather, 
we use the advantages of our scheme, particularly the two- 
dimensional nature of the code, to show that tidal forces 
(which are calculated self-consistently) can have a significant 
effect on the accretion flow in SXT discs and can account 
for the rebrightening observed in the X-ray light curves. 



6.2 Rebrightening time-scales 

The rebrightening in all the simulations begins 15-20 days 
after the onset of the outburst, while the peak of the second 
maximum occurs at different times according to the viscos- 
ity. This suggests that the rise-time of each rebrightening 
is controlled by viscous effects, but the time at which the 
rebrightening occurs is not. In each case, the accreted gas 
that is responsible for the rebrightening in the simulation 
originates from the outer (initially unirradiated) part of the 
disc. 

The viscous time-scale for the outer part of a Shakura- 
Sunyaev Q-disc is longer than 50 days. The viscous time- 
scale is 



tv 



aCsH' 



(27) 



100 




1.0 

Radius (10" cm) 



10.0 



Figure 10. Viscous time-scales. The solid lines are tvisc for an 
Q-disc with Q = 0.1 and Mi = 4M0. The upper solid line is 
for M = lO-iOMeyr-l, the middle M = IQ-^MQyr-l and the 
lower for M = 10~*MQyr~^. The dashed lines show tvisc in the 
SPH simulations, with Cs = 0.05aQb, Mi = 4M0, C = 0.5 (top), 
C = 1.0 (middle) and C = 1-5 (bottom). 



tv ~ 390 



Mf Rji d. 



(28) 



while the a-disc solution with a Kramer's opacity law gives 



where Mie is the mass accretion rate in units of 10 gs~ , Mi 
is the primary mass in solar masses and Rn is the radius in 
units of 10^^ cm. This is plotted in fig. |l^ for the parameters 
of A0620-003, along with the SPH viscous time-scale. 

Clearly, for the parameters of A0620-003 in outburst, 
the viscous time in the outer part of an Q-disc is longer 
than 50 days. We have seen in section 3 that for such an 
extreme mass ratio (q — 0.07), the 3:1 tidal resonance is 
easily accessible to the disc, and it is here that tidal effects 
start to affect the outburst. The tidal instability gives rise to 
superhumps observed in the V-band light curve of CVs, and 
we have previously shown for a CV with q = 0.15 that tidal 
effects prolong the outburst and drive eccentricity growth in 
the disc, with the launch of a two-armed spiral wave (Truss 
et al., 2001). Exactly the same situation applies to the tran- 
sient system studied here. Fig. 13 shows the evolution of the 
disc over run A - the spiral wave and the eccentricity are 
clearly evident during the outburst. 

The disc comes into contact with the 3:1 resonance soon 
after the onset of the outburst (flg. 13, top panel), but it 
takes some time for eccentricity to develop in the disc and 
the spiral waves to penetrate right down into the inner disc 
regions. Only then does the rebrightening occur (fig. 13, cen- 
tre panel). By this time the disc has developed a significant 
eccentricity and gas near the disc edge is spread over a range 
of different radii. So, what is the viscous time-scale near the 
disc edge? Consider the disc shown in the centre panel of fig. 
13. Gas at the edge of the disc nearest the secondary is at 
a radius R ~ 0.3a, while gas at the far edge is at R ~ 0.6a. 
This factor of two in radius translates to a factor of two in 
viscous time-scale. 

In a real disc there is an additional factor which could 
drive accretion on a shorter-than-expected time-scale and 
produce the rebrightening. Consider an idealised accretion 
disc with a Shakura-Sunyaev a model structure. The equa- 
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Figure 11. The analytical solution for mid-plane temperature 
in an accretion disc with Mi = S.TMq and Porb = 7.75h as 
a schematic demonstration of the action of tidal heating. The 
dashed line is calculated from equation ^ and does not include 
the tidal term. The solid line includes tidal torque calculated from 
equation ^ with n=5 and an arbitrarily large value of C. 



tion of energy balance in the unirradiated part of the disc, 
ba, is now: 



4crTn 



3kE 



9 ^GM ^ 



(29) 



where Qtido is the energy dissipation due to tides. This can 
be written in the form given by Papaloizou & Pringle (1977): 



Qtido — (^^Keplcr — f^orb)Ttidc 



(30) 



where Ttido is the tidal torque per unit area. In 1-D treat- 
ments of disc instabiliies, this is usually included as 



Ttidc — CRj/Snorb ( — 



(31) 



(Smak 1984; Hameury et al. 1998). There are two uncertain- 
ties in this equation: the value of C, which is a numerically- 
determined constant, and the value of n, which determines 
how much of the disc is tidally heated. Nevertheless, we can 
use this to demonstrate the effect of tides on the mid-plane 
temperature of the disc. This is shown schematically for n=5 
and an arbitrary value of C in fig. |l^. Clearly, in order to 
significantly impact on the mid-plane temperature, the tidal 
torque on the outer disc has to be both large and sustained. 

Figure |l^ shows the energy dissipation rate in the disc 
averaged over one orbital period near the peak of the simu- 
lated outburst in run A, along with a similar plot for a sim- 
ulated outburst of the dwarf nova SS Cygni (with the same 
sound speed) . SS Cygni is not tidally unstable, since the disc 
never encounters the 3:1 radius - the dissipation rate falls off 
with radius to the edge of the disc and rises in a well-defined 
narrow region associated with the stream-impact hot-spot. 
However, in A0620-003, tidal forces produce the dominant 
contribution to dissipation over a surprisingly large area of 
the disc. The tidal heating is spread over a range of radii 
because the disc is eccentric - over an orbital period the disc 
edge is not confined to a single radius. The dissipation rate 
in the outer regions of the disc is consequently about an or- 
der of magnitude higher than expected. We stress that this 
enhanced dissipation does not have any effect on the viscous 
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Figure 12. Normalised energy dissipation rate averaged over one 
orbital period in outburst simulations of A0620-003 (solid line) 
and the cataclysmic variable SS Cygni (dashed line). SS Cygni is 
not subject to the tidal instability - the increased dissipation in a 
narrow region near the disc edge is due to the stream-impact hot- 
spot. In stark contrast, the tidal heating in A0620-003 is totally 
dominant, with a local dissipation rate over an order of magnitude 
larger than the intrinsic disc dissipation rate acting over a very 
large area of the disc. The dissipation rates have been normalised 
to the rate at R = 0.15Rl(1). The 3:1 radius of A0620-003 is 
marked with a vertical line. 
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Figure 13. Evolution of surface density throught the outburst 
in run A. N is the total particle number in the disc. The curves 
correspond to t = 2 days (solid line), t = 30 days (dashed line) 
and t = 60 days (dot-dashed line) . Most of the accreted mass goes 
into producing the first maximum in the light curve (fig. the 
disc loses over 60% of its mass in the first 20 days. The secondary 
maximum is caused by the arrival of mass initially near the disc 
edge (the large peak in the solid curve). 



time-scale in the simulations : the viscosity is fixed at the 
value set in equation |l7| . 

Figures ^ and |14| show the evolution of surface den- 
sity and particle position in the disc. The disc edge very 
quickly encounters the 3:1 resonant radius and a two- armed 
spiral wave is launched (fig. |l4[ top panel). At the peak of 
the outburst, material accumulates in the spiral arms in the 
unirradiated outer regions of disc (solid curve in fig. |l3| ), 
while gas is the inner region is accreted. The surface density 
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Figure 14. Evolution of the disc through the outburst. The plots 
from top to bottom correspond to the same times as the three 
density profiles in figure hsi and are centred on the black hole 
primary. Particles are injected from the Li point to the left of the 
plots; the inner part of each disc (white) is not modelled. Near the 
peak of the primary maximum (top), Rjj = 0.35a and the edge 
of the disc has encountered the 3:1 resonance (R = 0.47a). The 
rebrightening is associated with the arrival of gas in the tidally- 
induced spiral arm which penetrates into the inner regions of the 
disc (centre). Successive waves of material are accreted in this 
way as the outburst progresses, but the variations in the X-ray 
light curve diminish as the disc is drained (bottom). 
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Figure 15. Accretion rate as a function of time for a simulation 
of a transient with the parameters of A0620-003, except that the 
mass ratio q = 0.6. The disc does not become tidally unstable, and 
there is no rebrightening, just a long tail. As before, fjjot = 1.0. 

in the outer 'spike' of gas eventually surpasses the upper 
critical limit to trigger the disc instability. As the outburst 
progresses, the disc eccentricity grows. The rebrightening is 
associated with the arrival of the initially unirradiated gas at 
the black hole. This is not viscously spread on a time-scale 
that would apparently be appropriate to the outer disc - 
the spiral wave launched by the tidal instability penetrates 
the inner part of the disc and carries with it some of the 
gas in the outer disc (centre plot in fig. Hence the re- 

brightening is caused by the accretion of a genuinely short 
pulse of gas. This penetration of the spiral arms into the 
inner regions of a hot disc has also been found in grid-based 
calculations performed by Stehle (1999). Successive waves 
of material are accreted in this way as the outburst pro- 
gresses. We find that ~ 60% of the disc mass is accreted 
in the primary maximum. This figure would be higher in a 
calculation in which a greater fraction of the disc is exposed 
to the radiation. 

The enhanced dissipation in the outer regions of the disc 
is consistent with the findings of Blondin (2000), who used 
an Eulerian scheme to show that tidal forces alone can pro- 
duce a viscosity of the order of a ~ 0.1 in the outer regions 
of an accretion disc. The effective viscosity decreased rapidly 
with radius, and only made a significant contribution in the 
outer half to two-thirds of the disc: exactly what we find 
with SPH in fig. |l^. It is still not absolutely clear that this 
will be sufficient to account for the 10 to 20 day time-scale 
for the rebrightening in A0620-003. Fig. |lo| shows that the 
viscous time-scale at the inner edge of the region of enhanced 
dissipation (R ~ 0.5Rl(1) ^ 6.5 x lO'^^cm) is ~ 85 days for a 
pure a-disc. However, Blondin (2000) and Spruit (1987) note 
that the propagation of spiral waves in a real disc is more 
likely than in an isothermal disc, because in a real disc the 
sound speed increases with decreasing radius (cs ^ R~s). 
Hence, with a more realistic equation of state, the effective 
Q may be even higher in the inner disc. Furthermore, there 
is an additional effect at work in the SPH simulations which 
was not found by Blondin - the development of a significant 
eccentricity in the disc. These two effects may combine to 
make the short time-scale of the rebrightening achievable. 

In the SPH simulations, both sound speed and viscos- 
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ity are fixed, so the rebrightening is caused by the arrival 
at the primary of a density wave associated with the spiral 
waves penetrating deep into the inner part of a highly ec- 
centic disc. We perform one further simulation to show that 
the rebrightening is indeed a tidal effect. The same input 
parameters as run A are used, with the exception that the 
mass ratio, q = 0.6. At this mass ratio, the 3:1 resonance 
should not be accessible to the disc and therefore it cannot 
become tidally unstable. The result is shown in fig. As 
expected, no significant eccentricity develops and there is no 
rebrightening: just a simple decay with a long tail. 

In summary, then, there are two time-scales associated 
with the rebrightening: the growth time-scale of the tidal 
instability in the disc, giving the time delay from the onset 
of the outburst to the start of the rebrightening, and the 
viscous time-scale of this tidally-stressed gas. In the simu- 
lations, the eccentricity of the disc causes gas initially near 
the disc edge to lie at much smaller radii. This explains why 
the rise-time of the rebrightening is so short : the gas that 
produces it is viscously spread on a time-scale much shorter 
than would be expected for gas originating near the edge of 
the disc. Undoubtedly, with a larger number of particles, the 
spatial resolution of the spiral shocks would improve and the 
simulated rebrightening would be even narrower. We have 
also proposed that in a real disc (but not the simulations) 
there is an additional factor that can increase the viscosity 
and make the rise-time of the rebrightening much shorter 
than expected. This is the tidal heating of the disc, raising 
the mid-plane temperature. 



The main limitation of this approach is that the viscous 
evolution of the disc is governed by the local surface den- 
sity and not a local mid-plane temperature. The ultimate 
goal is to merge the current ID thermodynamic state-of- 
the-art with our 2D approach. The eventual extension to 
three dimensions is also particularly desirable in examining 
the beahviour of the tidal forces and the dynamics of the 
spiral arms. 

There already exists a wealth of X-ray observations of 
transient sources, and we stress that although some aspects 
of the X-ray light curves are found to be common to many 
SXTs, they are all different to some degree. In this work 
we have presented a mechanism for the outbursts which ex- 
plains the origin of the rebrightenings observed in the X-ray 
light curves of some of these systems. 

Although one would not expect to observe superhumps 
in the X-ray light curve, superhumps have been seen in op- 
tical light curves of SXT outbursts (O'Donoghue & Charles 
1996). The superhump mechanism must be slightly differ- 
ent to that operating in cataclysmic variables, as the optical 
emission from the disc is completely dominated by the re- 
processed X-rays. Rather, the superhump modulation can 
be attributed to the change in visible area of an eccentric 
disc over an orbital period (Haswell et al. 2001). This is ex- 
actly the scenario which we find in this work, and is a clear 
demonstration of the importance of using a two- (or prefer- 
ably a three-) dimensional code able to take the tidal effects 
into account in a natural way. 



7 CONCLUSIONS 

We have simulated complete outbursts of soft X-ray tran- 
sient systems, for the first time in two dimensions. The re- 
brightenings seen in the X-ray light curves of several systems 
can be explained by competition between the ionization in- 
stability, the effects of disc irradiation and tidal forces from 
the secondary star. Their appearance in the light curve is 
ultimately determined by the irradiated area of the disc, 
the location within the disc where the disc instability is first 
triggered and the growth time-scale of the tidal instability. 
The main findings of this work can be summarized as: 

• A rebrightening can occur if a small portion of the disc 
is shielded from the central X-rays in some way. This com- 
plements the conclusions of Dubus et al. (1999), who showed 
that a planar disc cannot be kept heated by irradiation from 
the central source, in contradiction to observations. 

• Tidal torques can be a dominant factor in the energy 
balance over a large area of the disc. Such an enahnced dis- 
sipation in a disc with the classical a-model structure will 
drive accretion on a shorter than expected time scale. 

• The time delay between the onset of the outburst and 
the start of the rebrightening is determined by the growth 
time-scale of the eccentricity and spiral waves in the disc 
associated with the tidal instability. 

• The rise-time of the rebrightening is governed by the 
disc viscosity - but tidal eccentricity effects cause this time- 
scale to be much shorter than expected for the arrival of gas 
from the disc edge. 

• No rebrightening is reproduced if the entire disc is irra- 
diated. 
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